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Abstract 

The designation of the customary restricted three-body disturbing potential as the 
perturbation Hamiltonian is believed to be the cause of Neptune ring arcs’ radial offset 
between theories and observations. To identify the appropriate perturbation Hamiltonian, 
the energy integral in the fixed frame of a restricted three-body system, consisting of the 
central, primary, and test bodies, is reconsidered. It is shown that the perturbation energy 
includes the disturbing potential $ and the potential arising from the angular momentum 
terms of the test body. Both potentials happen to be singular as the test body goes to 
infinity contradicting to the perturbation nature. These two potentials can be combined 
to an energy relevant disturbing potential = ( 3 <f> which is regular at infinity because of 
the cancellation of the singularities. For circular orbits of the primary, the energy equa¬ 
tion is conservative, and <!>* is identified as the perturbation Hamiltonian. Applying this 
result to evaluate the backgrund effect of Triton to the arc-Galatea system of Neptune, 
it is shown that there is a small difference A<h = (<1>* — <h) which amounts to an outward 
radial offset of the corotation location of Galatea by 0.3 Km. The mismatch between 
the pattern speed of Galatea’s corotation potential and the mean motion velocity of the 
arcs could be resolved by considering the finite mass of Fraternite. However, by using $*, 
Galatea’s eccentricity could be reassessed in terms of the mass of Fraternite. 
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1 Introduction 

Since the discovery of the Neptune ring arcs [Hubbard et al 1986], the constant monitoring of 
their evolution have indicated that these arcs are dynamically stable [Smith et al 1989, Sicardy 
et al 1999, Dumas et al 1999]. The ring arcs are located at a radius r s from Neptune slightly 
larger than the radius tq of a nearby known satellite Galatea. Current theories to account for 
the arcs fall in two scenarios. The first is the multiple moons scenario that considers one or 
more unidentified Lagrange moons at r s to azimuthally confine the arcs at L4, L5 Lagrangian 
points of the arc-moon-Neptune restricted three-body system, and the shepherd moon Galatea 
at t'g to do the radial confinement [Lissauer 1985, Sicardy and Lissauer 1992], In contrast to 
the multimoon theories, the second is the single moon scenario that has Galatea as the only 
body to provide a stable configuration to the ring arcs through corotation-Lindblad 42/43 orbit- 
orbit resonance via the disturbing potential derived from the arc-Galatea-Neptune restricted 
three-body system [Goldreich et al 1986, Porco 1991, Foryta and Sicardy 1996]. The corotation 
potential is obtained by expanding the disturbing potential in terms of the orbit elements of 
Galatea, and by taking the time average to remove the fast oscillating non-resonant terms. The 
Lindblad potential, which describes the interaction of the test mass with the corotation site, is 
reached by expanding the disturbing potential in terms of the orbit elements of the test mass. 
With the aid of orbital measurements of the arcs, an analysis on the location of the corotation- 
Lindblad resonance was developed by Horanyi and Porco [1993]. The role of the massive 
satellite Triton on the 42/43 orbit-orbit resonance can be included by linearly superpositioning 
the arc-Triton-Neptune restricted three-body effects. Nevertheless, recent precise observations 
have shown a small but very critical orbit offset of sbout 0.3 Km from the predicted radial 
position [Sicardy et al 1999, Dumas et al 1999, Showalter 1999, Namouni and Porco 2002] that 
can not be resolved by fine tuning the system parameters. 

In order to resolve this corotation resonance orbit offset, we first summarize the energy 
equation of the restricted three-body system consisting of the central, primary, and test bodies 
in Sec.2. This energy integral in the fixed frame contains the kinetic and potential energies on 
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the left side and the energy exchange terms on the right side. The potential energy contains the 
central force potential of the central body, the disturbing potential of the primary body, and a 
potential arising from the angular momentum of the test body. The disturbing potential gives 
the disturbing force in the equations of motion through its negative gradient, and is customary 
identified as the perturbation Hamiltonian in analytic orbit-orbit interaction models. This 
designation is inappropriate, and we believe it is the cause of the radial offset for two reasons. 
First, the disturbing potential does not represent the entire perturbed energy. There is the 
potential energy corresponding to the angular momentum of the test body that is left out. 
Second, the disturbing potential as it is diverges at infinity although the disturbing force is 
regular there. It is also noted that the angular momentum potential of the test body diverges 
as well at infinity. 

We next seek to express the energy integral in an alternative form. In Sec.3, the angular 
momentum term on the left side and the energy exchange terms on the right side of the energy 
equation are evaluated to the zeroth order of eccentricity of the primary body. These con¬ 
tributions are shown to be proportional to the usual force relevant disturbing potential itself. 
Grouping the disturbing potential terms together defines an energy relevant disturbing poten¬ 
tial which puts the energy integral invariant to zeroth order of eccentricity. This energy relevant 
disturbing potential is a regular function because the singularity of the usual disturbing poten¬ 
tial is canceled by the singularity of the angular momentum contribution. The perturbation 
Hamiltonian based on the energy relevant disturbing potential is constructed in Sec.4. 

This energy related perturbation Hamiltonian is applied to the neptunian system to evaluate 
the arc offset. The background effect of Triton is superimposed on the arc-Galatea-Neptune 
three-body system. Since the energy related disturbing potential of Triton differs slightly from 
the force related one, a small contribution is added to the corotation potential of Triton which 
amounts to the observed outward radial offset of the arcs. The requirement on Galatea’s residue 
eccentricity and the mass of the main arc Fraternite to close the mean motion mismatch are 
estimated anew. 
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2 Energy Integral-The Jacobi Constant 

Here, we consider a restricted three-body system with a central body S, a primary body X, 
and a test body s where M, m x are the masses of the central body S and the primary body X 
respectively. Also r x = (■ r x , 9 X ), f s = (r s , 9 S ) are the position vectors of X and s measured from S 
with respect to a fixed reference axis in space. Furthermore, R = r s — r x is the position vectors 
of s measured from X. We consider all the bodies moving on the ecliptic plane by neglecting 
the orbit inclinations. With respect to a coordinate system centered at the central body S, the 
equations of motion are 
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Here, A 9 S = (9 S — 9 X ), whereas u s = d9 s /dt and u> x = d9 x /dt are the angular velocities of 
s and X about the central body S respectively with respect to a fixed reference axis. The 
right sides of Eq.(3) and Eq.(4) are the radial and transverse disturbing forces on s, and 
can be expressed in terms of the gradient of the disturbing potential given by v4 > (fb, r a; ) = 
(d$/dr s )r s + (l/r s )(d$/d6 s )6 s . 
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Because of the presence of the primary body, the force is time dependent. The energy 
difference between two points in space is dependent on the trajectory of s, and the force held 
is not conservative even though it is derivable from a potential <f>. This disturbing potential $ 
is often identified as the perturbation Hamiltonian in analytic models. Nevertheless, we should 
point out that the second term in <f>, the indirect term, is r s dependent and it diverges at 
infinity. However, this singularity in the potential does not affect numerical integrations of the 
equations of motion since the disturbing force is regular. It only concerns models that base on 
perturbation Hamiltonian where energy, not force, of the system is regarded as a fundamental 
functional. 

To derive the perturbation Hamiltonian, we consider the energy integral of s. We multiply 
Eq.(3) by dr s /dt and with standard manipulations, the energy equation reads [Tsui 2000, 2002] 


1 d ,dr s 2 

2 dt dt 


1 d , . o 

H--(r s o;J 

2 dt K J 


d r GM 1 Gm x 

—- + $ + -- 

dt r s 2 r x 


+ r 2 s u s u x 


= Hr 2 s"s) 


du x 

dt 


„ r/ 1 2 xr . . . . , 3 1 ,dr 

GW(— + yy)[r s cos(A0 s ) - r x \ + -—] 


2 r x 


dt 


= -(r&s) 


du x 

dt 


P. 


sX 


(5) 


On the left side of this equation are the kinetic and potential energies. The potential energy 
contains the contributions of the central body GM/r s , the disturbing potential <f> of the primary 
body, the angular momentum (r^ou s )u> x of the test body which is singular as r s becomes infinite. 
We have also taken the liberty to add the d{Gm x /2r x )/dt term on both sides of the equation 
so that the potential energy contains this added term. On the right side, there is an angular 
momentum term multiplied by the du x /dt. Since this term is multiplied over by du> x /dt , not 
u x , it is not a power exchange term of X, but a noninertial force term like centrifugal and 
Coriolis forces. Also, there is the second term, P s x, which is proportional to the radial velocity 
dr x /dt. This is a power exchange term of s with X. Due to the added term d(Gm x /2r x )/dt 
on both sides of the equation, the coefficient of this dr x /dt term vanishes at the Trojan points 
with R — r x — r s and A 6 S = tt/3. 



6 


In order to eliminate the time dependent nature of the disturbing potential, we transfer to 
the rotating frame of X. With respect to the SX reference axis, we write cu s = cu x + lv' s , where 
u' s is the relative angular velocity of the body s to the body X. We then have 


1 d ,dr s 2 

2 dt dt 


1 d 

2 dt' 


+ ~—{r s u s 
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d ,GM 


dt 1 


+ $ + - 

r, 2 
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+ 2 r s w i] 
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( 6 ) 


In this frame, X appears stationary and the angular momentum term becomes the also singular 
centrifugal r 2 cu 2 /2 term with the kinetic energy of s relative to the rotating X. The potential 
function now contains only spatial coordinates and is conservative. For the special case of 
circular orbits of X, the terms on the right side of Eq.(6) vanish. The integration constant 
of the now integrable Eq.(6) corresponds to the total energy of s which identifies $ plus the 
r 2 LU 2 /2 singular term as the perturbation Hamiltonian. 
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3 Disturbing Potential 

In orbit-orbit resonance interactions, representations in the fixed frame centered at S with a 
fixed reference axis in space is often preferred. We, therefore, reconsider Eq.(5) with special 
attention on the singular terms proportional to powers of r s . We combine the d(r^ou s uj x )/dt 
term on the left side with the (■ rgOU s )dou x /dt term on the right side of Eq.(5) to get 


1 d dr s 2 

2 dt dt 


Id. .9 

~t~- 7~ s^s) 
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This gives the first term on the right side which corresponds to the power provided by the 
torque on s applied by X. This is a power exchange term with X just like the P s x term. We 
now make an important remark that the angular momentum is an implicit function of the 
disturbing potential through Eq.(4), and it can be written in terms of to read uu x (d&/dO s ). 
Integrating the above equation over time gives 
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The first term on the right side corresponds to energy exchange of the motion in the tangential 
direction and the second term is in the radial direction. Since 6 S and 9 X appear in <f> only 
through A 9 S = (9 S — 9 X ), the derivative with 9 S in the first term is equivalent to differentiating 
with A 9 S . 

We note that these exchanges terms on the right side are evaluated along the trajectory of 
s. To evaluate these terms and to bring out the essential features, we expand them in powers 
of eccentricity of X. The zeroth power corresponds to circular orbit energy exchange with the 
second term on the right side P s x = 0 . As for the first term, with circular orbits of s and 
A 9 S = (ca s - Lu x )t, we get 
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The higher power terms in eccentricity constitute a residue energy exchange and could be ne¬ 
glected under circular orbits of X. Transporting this term to the left side and combining with 
the existing disturbing potential gives the Hamiltonian 
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This Hamiltonian of Eqs.(10) is valid for corotational s with X. ft is not valid for coorbital s 
with X because the relative motion of s is very slowly varying that A 6 S and r s in the disturbing 
potential are varying on the same slow time scale that the evaluation of the time integral on 
the right side of Eq.(8) is not justified. In this coorbital case, it is natural to use the rotating 
frame representation of Eq. (6) where the external torque on s is represented by the centrifugal 
potential. 

To summarize our results, we have reconsidered the angular momentum terms on both 
sides of Eq.(5), and evaluated their circular orbit contributions to the energy integral. We 
have recasted the energy integral in the form of Eqs.(10) where the perturbed energy is clearly 
identified by <h* = /3Q. This energy relevant disturbing potential <f>* is the sum of the force 
relevant usual disturbing potential $ in Eqs.(3) and (4) and the angular momentum (' r 2 s uj s ) 
terms of Eq.(5). Both potentials are singular as r s becomes infinite. It happens that the 
singularity of one cancels the singularity of the other making <f>* = /3<E> regular and well defined 
at infinity. Should s be well inside the orbit of X, (3 is slightly larger than unity. Should s be 
well outside the orbit of X, it is negative but almost null. In the case of s corotational with a 
nearby X, the factor f3 could enhance <f> several times increasing the strength of capture into 


resonance. 
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4 Corotation-Lindblad Resonance 

We now consider the Neptune arcs’ confinement according to the 42/43 Lindblad-corotation 
resonance with Galatea under the restricted three-body framework where an arc s as the test 
body, Galatea G as the primary body, and Neptune S as the central body. Since we are using 
a coplanar model, the corotation resonance is of eccentricity origin. Although this may not 
correspond to the actual orbit-orbit interaction of this system, this detail is not essential to 
our primary objective which is to account for the radial offset. By expanding the disturbing 
potential of s, Eqs.(10), to first order in eccentricity of G, with a G = a G /r < 1, and under the 
convergence conditions of Taylor expansion to first order in eccentricity 


e G a G 
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•C 1 , 2 ne G <C 1 
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the secular and resonant potentials of corotation (C) resonance are 
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where L% = GMr G , (3 G = [1 + uj g /(uj s - u G )], and 0 nGC out = [nO - (n - 1 )9 G - <j> G ] is the 
sG corotation resonance variable with s outside G. The subscript sGSC refers to the relevant 
three-body sGS system in minor, primary, central body order and the resonant mode C for 
corotation. A test mass would feel a time independent potential if its motion is such that 
<pnGCout is constant. Taking a time derivative on (f>n G Cout , therefore, would give 


d(f> n GCout r dd dO G d(j) G 
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This gives n spatially periodic potential sites set up by the primary body G for the test mass 
s when the mean motion velocity dd/dt of the test body fulfills the condition above. This 
velocity dd/dt is the pattern speed of the corotation resonance potential. Likewise, expanding 
the disturbing potential of s to first order in eccentricity of s gives the secular and resonant 
potentials of Lindblad (L) resonance 


1 m G 
+ 2 M 
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l m G GM 2 (o). s 

2 lU tU 
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(13) 


where 4> nG Lout = [('«- + 1)0 — n0 G — 4>\ is the sG Lindblad resonance variable with s outside G. 
This is the contribution to the resonance potential by the eccentricity of the test body. 

The effects of Triton X on this Gs corotation-Lindblad resonance can be linearly superim¬ 
posed by nesting s and G into the Triton-Neptune XS background system where both Galatea 
G and arc s appear as test bodies. The disturbing potentials of G and s in this XS background 
then read 
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Under this background, the contribution of Triton to the corotation secular potential of G can 
be obtained by expanding the disturbing potential of G, Eq.(14a), to first order in eccentricity 
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of G. With a! = ac/r x , it reads 
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Superpositioning the two contributions and droping the subscript secular give the total corota¬ 
tion secular potential of G 
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Next, we consider the Lindblad (L) secular potential of s by expanding Eq.(llb) to first order 
in eccentricity of s. With a = r/r x , the total Lindblad secular potential of s is 
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5 Radial Offset and Mean Motion Mismatch 

Under the restricted three-body framework and the superposition of Triton effects with the 
usual disturbing potential <f>, Eq.(16) reads 


total 


1 tyiq , GM 


) 2 a G b ( y 2 (a G ) + 


1 m x 

2 ~M 



^sGSC + &GXSC ■ 


(18) 


This secular term defines a negative total potential function —[GM/r s + ^c,total]- The intersect 
of this potential function with the Hamiltonian constant of s gives the radial position of s 
which defines the expected resonance location. However, the observed location differs from this 
expected one by 0.3 Km outwards. 

We believe the orbit offset is due to the fact that the disturbing potential $ being singular 
and that the also singular angular momentum (■ r 2 u s ) terms in the energy equation being ne¬ 
glected. Recoursing to the procedured in earlier Sections, our disturbing potential <f>* = /3<f> of 
Eq.(lOb) is regular at infinity and with 


®C,total - Pg^sGSC + PxG^GXSC ■ 


(19) 


The difference of energy between these two approaches is 

A^ C ,total — (Pg — 1 )*^sGSC + (PxG ~ 1 ) < ^GX 5 'C ~ (PxG ~ l^G.YSC , ( 20 ) 


where the ajustment of Triton dominates over that of Galatea. Taking this difference into con¬ 
siderations, the potential function would be — [GM/r s +^ G ,totai+A^ G ,total] which is slightly more 
negative. The intersect with the Hamiltonian constant would, therefore, give a slightly outward 
location. As an estimate, total is given by (l/2)(m x /M)(GM/L x ) 2 [u x /(u G - u x )]bfj 2 {a'). 
Taking < E >= —(1/2 ){GM/r G ) introduces an offset 
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( 21 ) 


Considering M — lx 10 26 Ah?, m x = 2 x 10 22 Kg, ( ra/r x ) = 1/6, ( oj x /ojg ) = 1/15, ^ 1/2 ( Cfc/ ) = 2, 
we get <5r/r G = 4.4 x 10 -6 . Taking r G = 2.5 x 25 x 10 3 Km gives Sr = 0.3 Km which is 
compatible to the observed offset [Namouni and Porco 2002], 

Another question with the Neptune arcs is the mean motion velocity mismatch of the arcs 
with the pattern speed of the corotation resonance potential of Galatea. According to recent 
measurements, the corotation resonance variable 4> n GCout = [nO — (n — 1 )6q — 4>g] is n °t exactly 
constant in time because of the mismatch. It has been suggested to consider the finite arc mass 
which would add to Galatea’s epicycle frequency [Namouni and Porco, 2002], Under a general 
three-body framework, with finite m s , the disturbing potential of s acting on G is 

CtTT) CtTD 
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K T'i 


Expanding to first order in eccentricity of G gives the Lindblad (L) resonance of G relative to s 
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where 4> n sLin = [(n — 1 )0g — nd + 0 G ] = — (j) n GCout is the Lindblad resonant variable of G with 
G inside s. The added contribution of s on the epicycle frequency of G can be calculated by 


2 d0 G 
a G u; G e G — 
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de G 
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To close the mismatch, Namouni and Porco [2002] have estimated that e G = 10 4 should 
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m s /mG = 0.26, or ea = 10~ 6 should m s /mG = 0.002. Nevertheless, ec could not be arbitrarily 
small, because it would eliminate the corotation resonance potential amplitude of eccentricity 
type. In such cases, we have to reach out for the 84/86 potential of inclination type whose 
angular width disagrees with the main arc, Fraternite, spread of 10 degrees. With onr energy 
related disturbing potential, we would have §* mis = Pmis^mis, and Eq.(24) would read 


2 d<pG 

a G LJG6G ^r 


d 

de G 


4 > 


* 

mis,resonant 
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d 

de G 
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mis,resonant 


(25) 


This would reassess the parameters. With /3 mis = 44, we could have = 10 4 should m s /rriG = 
0.006. 
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6 Conclusions 

We have reconsidered the disturbing potential <f> of the restricted three-body system. While 
the disturbing force in the equations of motion is regular, we have noted that the disturbing 
potential gets very large at large distances and becomes singular at infinity. This characteristics 
is in contradiction to the potential as a small perturbation. Together with this singularity, there 
are angular momentum (r^u s ) terms on both sides of the energy integral that are also singular 
at infinity. These angular momentum terms can be casted in terms of <f> itself which leads to 
an energy equation with <f>* = f3 <f>. Because of the mutual cancelation of singularities from <f> 
and (r^oj s ), the energy related disturbing potential <f>* is regular and converges at infinity. 

For circular orbits, the energy equation takes a conservative form, and the perturbation 
Hamiltonian can be identified by <£>*. Comparing to the perturbation Hamiltonian given by <f> 
from the non-conservative energy equation, there is a difference A$ = (<f>* — <£>) = (/? — 1)$. 
To understand this difference which removes the singularity at infinity, we note that the often 
neglected angular momentum terms are implicitly dependent on <f> through Eq.(4), although 
they are not explicit functions of <f>. 

As a test case of this regular and analytic disturbing potential <£>*, we have applied to the 
corotation-Lindblad resonance of Galatea-arcs configuration in the Neptune system with the 
background effect of Triton linearly superpositioned. It is shown that the radial offset of the arcs 
could be accounted for by the difference of the secular contributions of the Triton background 
A Triton- As for the mismatch between the pattern speed of Galatea’s corotation resonance 
potential and the mean motion velocity of the arcs, we have followed the approach of a massive 
arc that enhances the epicycle velocity of Galatea. Nevertheless, by using the regular analytic 
disturbing potential <f>* = /?$ for interactions, the eccentricity of Galatea could be reassessed. 



16 


Acknowledgments 

This work was supported by the Conselho Nacional de Desenvolvimentos Cientifico e Tec- 
nologico (CNPq, The Brazilian National Council of Scientific and Technologic Developments) 
and the Fundacao de Amparo a Pesquisa do Estado do Rio de Janeiro (FAPERJ, The Research 
Fostering Foundation of the State of Rio de Janeiro). 



17 


References 

1- Dumas, C., R.J. Terrilc, B.A. Smith, G. Schneider, and E.E. Berklin 1999. Stability of 

Neptune’s ring arcs in question. Nature 400, 733-735. 

2 - Foryta, D.W., and B. Sicardy 1996. The dynamics of the Neptunian Adams ring’s arcs. 

Icarus 123, 129-167. 

3- Goldreich, P., S. Tremaine, and N. Borderies 1986. Towards a theory for Neptune’s arc 

rings. Astron. J. 92, 490-494. 

4 - Horanyi, M., and C.C. Porco 1993. Where exactly are the arcs of Neptune? Icarus 106, 

525-535. 

5- Hubbard, W.B., A. Brahic, B. Sicardy, L.R. Elicer, F. Roques, and F. Vilas 1986. Occul- 

tation detection of a neptunian ring like arc. Nature 319, 636-640. 

6- Lissauer, J.J. 1985. Shepherding model for Neptune’s arc ring. Nature 318, 544-545. 

7 - Namouni, F., and C.C. Porco 2002. The confinement of Neptune’s ring arcs by the moon 

Galatea. Nature 417, 45-47. 

8- Porco, C.C. 1991. An explanation for Neptune’s ring arcs. Science 253, 995-1001. 

9 - Showalter, M.R. 1999. Neptune’s misbehaving rings. Nature 400, 709-710. 

10 - Sicardy, B., and J.J. Lissauer 1992. Dynamic model of the arcs in Neptune’s 63 K ring. 

Adv. Space Res. 12, 81-95. 

11 - Sicardy, B., F. Roddier, C. Roddier, E. Perozzi, J.E. Graves, O. Guyon, and M.J. Northcott 

1999. Images of Neptune’s ring arcs obtained by a ground based telescope. Nature 400, 
731-732. 

12 - Smith, B.A., et al 1989. Voyager 2 at Neptune: imaging science results. Science 246, 


1422-1449. 



18 


13 - Tsui, K.H. 2000. Asteroid-Planet-Sun interaction in the presence of a planetary satellite. 

Icarus 148, 139-146. 

14 - Tsui, K.H. 2002. Satellite capture in a four-body system. Planet. Space Sci. 50, 269-276. 



